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Abstract 

We present a measurement of the deuteron spin-dependent structure function gf based on 
the data collected by the COMPASS experiment at CERN during the years 2002-2004. 
The data provide an accurate evaluation for Tf, the first moment of gi{x), and for the ma- 
trix element of the singlet axial current, oq. The results of QCD fits in the next to leading 
order (NLO) on all gi deep inelastic scattering data are also presented. They provide two 
solutions with the gluon spin distribution function AG positive or negative, which describe 
the data equally well. In both cases, at = 3(GeV/c)^ the first moment of AG(x) is 
found to be of the order of 0.2 - 0.3 in absolute value. 
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The spin structure function gf of the deuteron has been measured for the first time almost 15 
years ago by the SMC experiment at CERN [1]. Since then, high accuracy measurements of gf 
in the deep inelastic scattering (DIS) region have been performed at SLAC [2,3] and DESY [4]. 
Due to the relatively low incident energy, the DIS events collected in those experiments cover 
only a limited range of x for > l(GeV/c)^, x > 0.015 and x > 0.03, respectively. Further 
measurements covering the low x region were also performed at CERN (see [5] and references 
therein). Besides its general interest for the understanding of the spin structure of the nucleon, 
gf is specially important because its first moment is directly related to the matrix element of the 
singlet axial vector current oq. A precise measurement of gf can thus provide an evaluation of 
the fraction of nucleon spin carried by quarks, on the condition that the covered range extends 
far enough to low x to provide a reliable value of the first moment. 

Here we present new results from the COMR\SS experiment at CERN on the deuteron 
spin asymmetry Af and the spin-dependent structure function gf covering the range 1 (GeV/c)^ < 

< 100(GeV/c)^ in the photon virtuality and 0.004 < a; < 0.7 in the Bjorken scaling vari- 
able. The data sample used in the present analysis was collected during the years 2002-2004 
and corresponds to an integrated luminosity of about 2 fb^. Partial results based on the data 
collected during the first two years of the data taking have been published in Ref . [6] . At the 
time, the values of gf were not precise enough, in particular at large x, to allow a meaningful 
evaluation of the first moment, Tf. The results presented here are based on a 2.5 times larger 
statistics and supersede those of Ref. [6]. We refer the reader to this reference for the description 
of the 160 GeV muon beam, the ^LiD polarised target and the COMR\SS spectrometer which 
remained basically unchanged in 2004. A global fit to all g^'"'''^ data is needed to evolve the 
gf{xi, Ql) measurements to a common Q"^. As previous fits were found to be in disagreement 
with our data at low x, we have performed a new QCD fit at NLO. The resulting polarised 
parton distribution functions (PDF) are also presented in this paper and discussed in relation 
with the new data, however without a full investigation of the theoretical uncertainties due, for 
instance, to the values of the factorisation and renormalisation scales. 

The COMPASS data acquisition system is triggered by coincidence signals in hodoscopes, 
defining the direction of the scattered muon behind the spectrometer magnets, and by signals 
in the hadron calorimeters [7]. Triggers due to halo muons are eliminated by veto counters in- 
stalled upstream of the target. Inclusive triggers, based on muon detection only, cover the full 
range of x and are dominant in the medium {x, Q^) region. Semi-inclusive triggers, based on 
the muon energy loss and the presence of a hadron signal in the calorimeters, contribute mainly 
at low X and low Q^. Purely calorimetric triggers, based on the energy deposit in the hadron 
calorimeter without any condition on the scattered muon, account for most events at large Q"^. 
The relative contributions of these three trigger types are shown in Fig. [T] as a function of x. The 
minimum hadron energy deposit required for the purely calorimetric trigger has been reduced 
to 10 GeV for the events collected in 2004. As a consequence, the contribution of this trigger 
now reaches 40% at large x, compared to 20% in 2002-2003 (Ref. [6]). 

All events used in the present analysis require the presence of reconstructed beam muon 
and scattered muon trajectories defining an interaction point, which is located inside one of 
the target cells. The momentum of the incoming muon, measured in the beam spectrometer, is 
centered around 160 GeV/c with an RMS of 8 GeV/c for the Gaussian core. In the present 
analysis its value is required to be between 140 and 180 GeV/c. In addition the extrapolated 
beam muon trajectory is required to cross entirely both target cells in order to equalize the fluxes 
seen by each of them. The scattered muon is identified by signals collected behind the hadron 
absorbers and (except for the purely calorimetric trigger) its trajectory must be consistent with 
the hodoscope signals defining the event trigger. For hadronic triggers, a second outgoing recon- 
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structed track is required at the interaction point. The DIS events used in the present analysis are 
selected by cuts on the four-momentum transfer squared {Q^ > l(GeV/c)^) and the fractional 
energy of the virtual photon (0.1 < y < 0.9). The resulting sample consists of 89 x 10^ events, 
out of which about 10% were obtained in 2002, 30% in 2003 and 60% in 2004. In order to ex- 
tend the coverage of the low x region, we also analyse events in the interval 0.003 < x < 0.004 
selected in the same way but with a cut lowered to 0.7 (GeV/c)^. These events are included 
in the figures but not used in QCD calculations or moment estimation, in view of their low Q^. 

During data taking the two target cells are polarised in opposite directions, so that the 
deuteron spins are parallel (H) or antiparallel (tJ,) to the spins of the incoming muons. The 
spins are inverted every 8 hours by a rotation of the target magnetic field. The average beam 
and target polarisations are about —0.80 (—0.76 in 2002 and 2003) and ±0.50, respectively. 

The cross-section asymmetry = {a^^ — cr^^)/(a^-^ + cr^^), for antiparallel (t|) and 
parallel (ft) spins of the incoming muon and the target deuteron can be obtained from the 
numbers of events Ni collected from each cell before and after reversal of the target spins: 

Ni = ai(Pin,W{l + PBPTfA''), z = 1,2,3,4, (1) 

where is the acceptance, (pi the incoming flux, nj the number of target nucleons, a the spin- 
averaged cross-section, Pb and Pt the beam and target polarisations and / the target dilution 
factor. The latter includes a corrective factor p = cr^i^ /cr^f^ [8] accounting for radiative events on 
the unpolarised deuteron and a correction for the relative polarisation of deuterons bound in ^Li 
compared to free deuterons. Fluxes and acceptances cancel out in the asymmetry calculation on 
the condition that the ratio of the acceptances of the two cells is the same before and after spin 
reversal [9]. 

The longitudinal virtual-photon deuteron asymmetry, Af, is defined via the asymmetry of 
absorption cross-sections of transversely polarised photons as 

Ai = {a^-aJ)/{2a^), (2) 

where aj is the 7*-deuteron absorption cross-section for a total spin projection J and is the 
total transverse photoabsorption cross-section. The relation between Af and the experimentally 
measured A'^ is 

A'' = D{Ai + r]Ai), (3) 

where D and r] depend on kinematics. The transverse asymmetry A2 has been measured at 
SLAC and found to be small [10]. In view of this, in our analysis, Eq. ^ has been reduced to 
Af ~ A'^/D. The virtual-photon depolarisation factor D depends on the ratio of longitudinal 
and transverse photoabsorption cross sections R = /a^. In the present analysis an updated 
parametrisation of R taking into account all existing measurements is used [11]. The tensor- 
polarised structure function of the deuteron has been measured by HERMES [12] and its effect 
on the measurement of the longitudinal spin structure was found to be negligible, which justifies 
the use of Eqs ([l]-[3]) in the present analysis. 

In order to minimize the statistical error of the asymmetry, the kinematic factors /, D and 
the beam polarisation Pb are calculated event-by-event and used to weight events. A parametri- 
sation of Pb as a function of the beam momentum is used, while for Pt an average value is 
used for the data sample taken between two consecutive target spin reversals. The obtained 
asymmetry is corrected for spin-dependent radiative effects according to Ref. [13]. The asym- 
metry is evaluated separately for inclusive and for hadronic events because the dilution factors 
and the radiative corrections to the asymmetry are different. This is because the correction due 
to radiative elastic and quasi-elastic scattering events only affects the inclusive sample. 
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It has been checked that the use of hadronic triggers does not bias the inclusive asymme- 
tries. The most critical case is for the calorimetric trigger events at large x, where high-energy 
hadron production is limited by kinematics. This effect has been studied by Monte Carlo, using 
the program POLDIS [14]. DIS events were generated within the acceptance of the calorimetric 
trigger and their asymmetry calculated analytically at the leading order. A selection based on 
the hadron requirements corresponding to the trigger was applied and the asymmetries for the 
selected sample compared to the original ones. The differences were found to be smaller than 
0.001 in all intervals of x (Fig.[2l) and thus negligible, so that inclusive and hadronic asymme- 
tries can be safely combined for further analysis (see also the SMC analysis [5]). 

The final values of Af(x, Q"^), obtained as weighted averages of the asymmetries in the 
inclusive and hadronic data sets, are listed in Table [H with the corresponding average values of 
X and Q^. They are also shown as a function of x in Fig. [3] in comparison with previous results 
from experiments at CERN [5], DESY [4] and SLAC [2, 3]. The values of Af confirm, with 
increased statistical precision, the observation made in Ref. [6] that the asymmetry is consistent 
with zero for x < 0.03. Values of Af originating from experiments at different energies tend to 
coincide due to the very small dependence of Af at fixed x. 

The systematic error of Af contains multiplicative factors resulting from uncertainties on 
Pb and P^, on the dilution factor / and on the ratio R = a^/a'^ used to calculate the depolari- 
sation factor D. When combined in quadrature, these errors result in a global scale uncertainty 
of 10% (Table[2l). The other important contribution to the systematic error is due to false asym- 
metries which could be generated by instabilities in some components of the spectrometer. In 
order to minimize their effect, the values of Af in each interval of x have been calculated for 
184 subsamples, each of them covering a short period of running time and, therefore, ensuring 
similar detector operating conditions. An upper limit of the effect of detector instabilities has 
been evaluated by a statistical approach. The dispersion of the values of Af around their mean 
agrees with the statistical error. There is thus no evidence for any broadening due to time depen- 
dent effects. Allowing the dispersion of Af to vary within its two standard deviations we obtain 
an upper limit for the systematic error of Af in terms of its statistical precision: agyst < OAagtat- 
This estimation accounts for the time variation effects of spectrometer components. 

Several other searches for false asymmetries were performed. Data from the two tar- 
get cells were combined in different ways in order to eliminate the physical asymmetry. Data 
obtained with different settings of the microwave frequencies, used to polarise the target by 
dynamic nuclear polarisation, were compared. No evidence was found for any significant appa- 
ratus induced asymmetry. 

The longitudinal spin structure function is obtained as 

where F2 is the spin-independent deuteron structure function. The values of gf listed in the last 
column of Table [T] have been calculated with the F2 parametrisation of Ref. [5], which covers 
the range of our data, and the new parametrisation of R already used in the depolarisation 
factor. The systematic errors on gf are obtained in the same way as for Af, with an additional 
contribution from the uncertainty on F^. The values of x, gf{x) for the COMR\SS data and, 
for comparison, the SMC results [5] moved to the of the corresponding COMR^SS point 
are shown in Fig. HI The two curves on the figure represent the results of two QCD fits at NLO, 
described below, at the measured Q"^ of each data point. 

The evaluation of the first moment Tf{Q'^) = gf{x, Q'^)dx requires the evolution of all 
gi measurements to a common Q^. This is done by using a fitted parametrisation gl^^ix, Q"^), 
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so that 



9iix, Ql) = gi{x, Q'^) + fiff Ql) - g('\x, Q 



(5) 



We have used several fits of gi from the Durham data base [15]: Bliimlein-Bottcher [16], GRSV 
[17] and LSS05 [18], and we have chosen Ql = 3(GeV/c)^ as reference because it is close 
to the average of the COMPASS DIS data. The three parametrisations are quite similar in the 
range of the COMPASS data and have been averaged. The resulting values of g^ = {gi + g^) / 2 
are shown as open squares in Fig. [51 For clarity we now use g^ instead of gf because the 
correction for the D-wave state of the deuteron has been applied: 

g^ix,Q') = g'lix,Q^)/{l-1.5iUn) (6) 

with ud = 0.05 ± 0.01 [19]. It can also be seen in Fig.[5]that the curve representing the average 
of the three fits does not reproduce the trend of our data for x < 0.02 and therefore cannot be 
used to estimate the unmeasured part of g^ at low x. 

In view of this, we have performed a new NLO QCD fit of all gi data at > 1 (GeV/c)^ 
from proton, deuteron and '^He targets, including the COMPASS data. The deuteron data are 
from Refs. [2-5], the proton data from Refs. [2,4,5,20,21] and the ^He data from Refs. [22-25]. 

In order to optimise the use of the COMPASS data in this fit, all x bins of Table [H except 
the last one, have been subdivided into three intervals (Fig. |6l). The number of COMPASS 
data points used in the fit is thus 43, out of a total of 230. 

The fit is performed in the MS renormalisation and factorisation scheme and requires 
parametrisations of the quark singlet spin distribution AS(a;), non-singlet distributions Aq3{x), 
Aqs{x) and the gluon spin distribution AG{x). These distributions are given as an input at a 
reference Q'^ (= Ql) which is set to 3 (GeV/c)^ and evolved according to the DGLAP equations. 
The resulting values of gi (x, Q"^) are calculated for the (xj, Qf) of each data point and compared 
to the experimental values. 

The input parametrisations are written as 

x"''(l-x)^'-(l + 7fcx) 

= ^^7i 77' ' 

Jq (1 - xY^ (1 + ^kx)dx 

where AFk represents each of the polarised parton distribution functions AS, Ags, Agg and 
AG, and rfk is the integral of AF^. The moments, rf^, of the non-singlet distributions Ags and 
Agg are fixed by the baryon decay constants (F+D) and (3F— D) respectively, assuming SU(3)/ 
flavour symmetry. The linear term 7a; is used only for the singlet distribution, in which case the 
exponent I3g is fixed because it is poorly constrained by the data. This leaves 10 parameters 
in the input distributions. In addition, the normalisation of E155 proton data is allowed to vary 
within the limits quoted by the authors of Ref. [21]. 

The optimal values of the parameters are obtained by minimizing the sum 



N=^^y:'M^)-gT^{x,,Qi 



i=l 



(^{xi, Ql 



(8) 



Here the errors o are the statistical ones for all data sets, except for the proton data of El 55 
where the uncorrected part of the systematic error on each point is added in quadrature to the 
statistical one. In order to keep the parameters in their physical range, the polarised strange sea 
distribution As(x)+ As (x) = (l/3)(AS(x)— Agg(x)) is calculated at every step and required to 
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satisfy the positivity condition | As(x) \< s{x) at all Q'^ values. A similar condition is imposed 
on the gluon spin distribution AG{x). The unpolarised distributions s{x) and G{x) used in this 
test are taken from the MRST parametrisation [26]. This procedure leads to asymmetric errors 
on the parameters when the fitted value is close to the allowed limit. 

The fits have been performed with two different programs: the first one uses the DGLAP 
evolution equations for the spin structure functions [27], the other one, referred to in [28], uses 
the evolution of moments. The fitted PDF parameters are compatible within one standard devia- 
tion and the two programs give the same -probabilities. In each program the minimisation 
converges to two different solutions, depending on the sign of the initial value of the gluon first 
moment r]G'- one solution with AG > 0, the other one with AG < (Fig.[5l). The fitted distri- 
butions of gf{x) differ at low x but are both compatible with the data. The two additional data 
points at a; < 0.004 and > 0.7(GeV/c)^, not used in the fit, have too large statistical errors 
to provide a discrimination between the two solutions. The values of the parameters obtained in 
the fits with positive and negative AG are listed in Table [3] with their statistical errors and will 
be discussed below. 

The integral of gi in the measured region is obtained from the experimental values 
evolved to a fixed and averaged over the two fits. Taking into account the contributions 
from the fits in the unmeasured regions at low and high x we obtain (Table H]): 

rf (g2 = 3(GeV/c)2) = 0.050 ± 0.003 (stat.) ± 0.003 (evol.) ± 0.005 (syst.). (9) 

The second error accounts for the difference in Q"^ evolution between the two fits. The system- 
atic error is the dominant one and mainly corresponds to the 10% scale uncertainty resulting 
from the errors on the beam and target polarisations and on the dilution factor. 
For comparison, the SMC result [5] was 

SMC {Q^ = 10(GeV/c)2) = 0.021 ± 0.007 (stat.) ± 0.014 (evol.) ± 0.003 (syst.). (10) 

while our resuk evolved to = 10(GeV/c)2 is 0.051 ± 0.003 (stat.) ± 0.003 (evol.) ± 
0.005 (syst.). The difference between these two results reflects the fact that the COMPASS 
data do not support the fast decrease of gf (x, Ql = 3(GeV/c)^) at low x which was assumed 
in the SMC analysis, and thus force the fit to be different. In the COMPASS analysis, the part 
of obtained from the measured region represents 98% of the total value. This correction of 
only 2% has to be compared to a correction of about 50% with respect to the measured value in 
case of the SMC analysis [5]. 

is of special interest because it gives access to the matrix element of the singlet axial 
current Oq which, except for a possible gluon contribution, measures the quark spin contribution 
to the nucleon spin. At NLO, the relation between and reduces to 

rf (Q^) = ^ (l - + 0{a^)) (ao(Q^) + ^a^) . (11) 

From the COMPASS resuk on rf (Eq. (9)) and taking the value of ag from hyperon /3 
decay, assuming SU{3)f flavour symmetry (ag = 0.585 ± 0.025 [29]), one obtains with the 
value of as evolved from the PDG value as{m\) = 0.1187 ± 0.005 and assuming three active 
quark flavours: 

ao = 3(GeV/c)2) = 0.35 ± 0.03 (stat.) ± 0.05 (syst.). (12) 

The quoted systematic error accounts for the error from the evolution and for the experimental 
systematic error, combined in quadrature. 
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The relation between and can also be rewritten in order to extract the value of the 
matrix element ao in the limit Q"^ ^ oo. Here we will follow a notation of Ref. [30] introducing 
a "hat" for the coefficient Cf and ag at this limit: 

9 6b 

The coefficients and C^^ have been calculated in perturbative QCD up to the third order in 
as{Q') [30]: 

(g2) = 1 _ 0.33333(^) - 0.54959(^)' - 4.44725(^)' 

Cf^(Q2) = 3.5833(^)' - 20.2153(^)'. 

With as evolved at the same order, one obtains 

ao = 0.33 ± 0.03 (stat.) ± 0.05 (syst). (13) 

It should be noted here that the data have been evolved to a common on the basis of a fit 
at NLO only. However, the choice of a value close to the average of the data is expected to 
minimise the effect of the evolution on the result quoted above. Combining this value with as, 
the first moment of the strange quark spin distribution in the limit Q"^ ^ oo is, found to be 

(As + As)q2^oo = -(ao - CLs) = -0.08 ± 0.01 (stat.) ± 0.02 (syst.). (14) 

3 

As stated before, this result relies on SU (3) / flavour symmetry. A 20% symmetry breaking, 
which is considered as a maximum [29], would shift the value of As + As by ± 0.04. 

Previous fits of gi, not including the COMPASS data, found a positive AG'(x) and a 
fitted function gf{x) becoming negative for x < 0.025 at = 3(GeV/c)^, as shown by the 
dotted line in Fig. [5l The new COMPASS data do not reveal any evidence for a decrease of 
the structure function at limit x — 0. For our fit the data are still compatible with a positive 
AG, as shown by the full line in Fig. \5\ However in this case a dip at x ~ 0.25 appears in 
the shape of gf{x) for — > l(GeV/c)^. Its origin is related to the shape of the fitted AG{x), 
shown in Fig.|7](left). Indeed, the gluon spin distribution must be close to zero at low x, to avoid 
pushing gf down to negative values, and is also strongly limited at higher x by the positivity 
constraint \ AG{x) \ < G{x). The whole distribution is thus squeezed in a narrow interval around 
the maximum at a; ~ ac/iaG + Pg) — 0.25. 

In contrast, the fit with negative AG reproduces very well the COMPASS low x data with 
a much smoother distribution of AG{x) (dashed line on Fig. [5]) and without approaching the 
positivity limit (Fig. U\ right). The (1 + 'jx) factor in the singlet quark distribution is not used in 
this case because it does not improve the confidence level of the fit. 

Comparing the fitted parameters for AG positive and negative (Table [3]), we observe that 
the parameters of the non-singlet distributions Aq^^x) and Aqs{x) are practically identical. The 
value of r^s is slightly larger in the fit with AG < 0, as could be expected since in this case 
AS(x) remains positive over the full range of x: 

ri^ (q2 = 3(GeV/c)2) = 0.27 ± 0.01 (stat.) (AG > 0), (15) 

= 3(GeV/c)2) = 0.32 ± 0.01 (stat.) (AG < 0). (16) 
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We remind that in MS scheme rj^ is identical to the matrix element ao. 
The singlet moment derived from the fits to all gi data is thus: 

rys (g2 = 3(GeV/c)2) = 0.30 ± 0.01 (stat.) ± 0.02 (evoL). (17) 

Here we have taken the difference between the fits as an estimate of the systematic error and 
do not further investigate other contributions related to the choice of the QCD scale or the 
PDF parametrisations. The singlet moment obtained with COMPASS data alone (Eq. (12)) is 
slightly above this value and its statistical error is larger by a factor of 3. As stated before, the 
main uncertainty on the COMPASS result is due to the 10% normalisation uncertainty from the 
beam and target polarisations and from the dilution factor. The fact that the COMPASS data 
are on average slightly above the world average can already be detected by a comparison of the 
measured gf values to the curves fitted to the world data (Fig. [5]). Hence derived from the 
COMPASS value of Ff is found to be slightly larger than r/s. 

The polarised strange quark distributions, obtained from the difference between AS(x) 
and Aqs{x) are shown in Fig. [8l They are negative and concentrated in the highest x region, 
compatible with the constraint |As(x)| < s{x). This condition is indeed essential in the deter- 
mination of the Agg parameters which otherwise would be poorly constrained. 

Although the gluon distributions strongly differ in the two fits, the fitted values of their 
first moments are both small and about equal in absolute value \t]g\ ~ 0.2 — 0.3. We have 
also checked the stability of these results with respect to a change in as{m?z): when as{rn\) 
is varied by ±0.005 the values of rfc are not changed by more than half a standard deviation. 
In Fig. |9] the existing direct measurements of AG/G [31-33] are shown with the distributions 
of AG{x) /G{x) derived from our fits with G{x) taken from Ref. [26]. The HERMES value 
is positive and 2a away from zero. The measured SMC point is too unprecise to discriminate 
between positive or negative AG. The published COMPASS point, which has been obtained 
from a partial data sample corresponding to about 40% of the present statistics, is almost on the 
AG > curve but is only 1.3cr away from the AG < one, so that no preference for any of the 
curves can be given so far. It should also be noted that the measured values of AG/G have all 
been obtained in leading order QCD analyses. 

In summary, we have measured the deuteron spin asymmetry Af and its longitudinal spin- 
dependent structure function gf with improved precision at Q"^ > l(GeV/c)^ over the range 
0.004 < X < 0.70. The gf values are consistent with zero for x < 0.03. The measured values 
have been evolved to a common by a new fit of the world gi data, and the first moment has 
been evaluated at = 3(GeV/c)^ with a statistical error smaller than 0.003. From we have 
derived the matrix element of the singlet axial current ciq in the limit oo. With COMPASS 

data alone, at the order a^, it has been found that Oq = 0.33 ± 0.03 (stat.) ± 0.05 (syst.) and 
the first moment of the strange quark distribution (As + As)q2^^ = —0.08 ± 0.01 (stat.) ± 
0.02 (syst.). We also observe that the fit of world gi data at NLO yields two solutions with either 
AG(x) > or AG (a;) < 0, which equally well describe the present data. In both cases, the 
first moment of AG(a;) is of the order of 0.2-0.3 in absolute value at = 3(GeV/c)^ but the 
shapes of the distributions are very different. 
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X range 


(x) 




± 


gf 






[(GeV/(.:)-] 






0030-0 0035 


0033 


0.78 


003 ± 009 ± 004 


090 ± 240 ± 107 


0035-0 0040 


0038 


0.83 


-0 004 ± 007 ± 003 


-0 097 ± 183 ± 082 


A AA/1 A AAC 

U.UU4— U.UUj 


A AA/I A 

U.UU40 


1 1 r\ 
1 . iU 


U.UU4 ± U.UUy ± U.UU4 


n nfio -1- n 01 n -1- n nsn 
U.Uoz ± U.ziU ± u.uoy 


U.UUj— U.UUo 


A AACC 

U.UUj J 


1 00 


U.UUo ± U.UU / ± U.UUo 


n nfio -1- n 1 ap. -I- n nfio 
U.Uoz ± U.i40 ± U.UOZ 


A AA/: A AAQ 

U.UUo— U.UUo 


A Am A 

U.UU /U 


1 

L.jy 


n nno -u n nnc; -u n nno 
— U.UUz ± U.UUo ± U.UUz 


n nQ/1 _i- n nsf; -i_ n nQfi 
— U.Uc)4 ± U.UoO ± U.UoO 


U.UUo— U.UiU 


u.uuyu 


i.Oi 


n ni n 4- n nnfi -i- n nnQ 

— U.UiU ± U.UUD ± U.UUo 


n 1 Qo 4- n n7Q -i- n nQc; 
— u.ioy ± u.U/o ± u.Uoo 


U.U lU—U.UZU 


n ni /1 1 

U.U141 


^ 

Z. i J 


n nno -i- n nn/i -u n nno 
U.UUz ± U.UU4 it U.UUz 


n ni7-un nQQ -un ni/i 
U.Ui( ± U.Uoo ± U.Ui4 


A AOA A AOA 

U.UzU— U.UjU 


A r\OAA 
U.UZ44 


j.lo 


n nnQ -i- n nnfi -i- n nnQ 
U.UUo ± U.UUo ± U.UUo 


n ni7-i-n nQc; -un nic; 
U.Ui/ ± U.UoD ± U.UiD 


A AO A A A/1 A 

U.UjU— U.U4U 


A AO/I A 

U.Uj4d 


A Of^ 


n nnn -i- n nns -i- n nn/i 
U.UUy ± U.UUo ± U.UU4 


nn/ii 1 nnQf^. 1 nnifi 
U.U4i ± U.UoD ± U.Uio 


A A/1 A A A<^r\ 

U.U4U— U.UoU 


A A/1 QT 

U.U4o / 


j.oU 


n ni 7 J- n nns j- n nn/i 
U.Ui 1 ± U.UUo ± U.UU4 


n nc;/i -i- n nof; 1 n m 
U.U04 ± U.UzO ± U.UiZ 


0.060-0.100 


0.0765 


8.53 


0.058 ± 0.009 ± 0.007 


0.121 ±0.019 ±0.014 


0.100-0.150 


0.121 


12.6 


0.095 ±0.013 ±0.011 


0.123 ±0.017 ±0.014 


0.150-0.200 


0.171 


17.2 


0.123 ±0.020 ±0.014 


0.103 ±0.016 ±0.012 


0.200-0.250 


0.222 


21.8 


0.183 ±0.028 ±0.021 


0.106 ±0.016 ±0.012 


0.250-0.350 


0.290 


28.3 


0.216 ±0.030 ±0.024 


0.077 ±0.011 ±0.009 


0.350-0.500 


0.405 


39.7 


0.343 ± 0.049 ± 0.038 


0.055 ± 0.008 ± 0.006 


0.500-0.700 


0.566 


55.3 


0.626 ±0.112 ±0.075 


0.027 ± 0.005 ± 0.003 



Table 1 : Values of A\ and g'{ with their statistical and systematical errors as a function of x with the 
corresponding average values of x and Q^. The minimum cut is 1 (GeV/c)^ except for the first 
two points where it is lowered to 0.7 (GeV/c)^. These two data points are shown on the figures as 
complementary information but were not used in the fits. 



Multiplicative 
variables 
error, AA^"'* 


Beam polarization 




5% 


Target polarization 


(IPt/Pt 


5% 


Depolarization factor 


dD{R)/D{R) 


2-3 % 


Dilution factor 


df/f 


6 % 


Total 




^^mult ^ 0.1 


Additive 
variables 
error, I\Af'^ 


Transverse asymmetry 


v/p ■ AA2 


10-^ - 5 ■ 10-3 


Radiative corrections 




10-^ - 10-3 


False asymmetry 


Afalse 


< 0.4 • AAf 



Table 2: Decomposition of the systematic error of Ai into multipUcative and additive variables contri- 
butions. 
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AG > 




Prog. Ref. [27] 


Prog. Ref. [28] 




0.270 ± 0.014 


0.284 1 




I", ono + 0.074 

— U.OUC) _ Q Qyg 


rv nn/- + 103 

-0.226 1 ^ofoi 


/3s 


r, /%ri+ 24 

3.601 0:22 


n + 30 

3.69 1 ^oil 


7s 


-16.01}:^ 


-15.8 l^i 


Vg 


n oof^+ 0.049 


n 900 + 0.040 
U.zoo _ 0.053 


«G 


0.44 


011+ 0.42 
0.53 




10 (fixed) 


10 (fixed) 




-0.226 ±0.027 


-0 226 + 

U.^^U _ 0.027 




^•^"J- 0.10 


9 qo + 0.11 
^•'-'O - 0.10 


as 


3S+ 

U.OO_ 0.44 


45 + 

U.^O _ 0.43 




q qc+ 0.60 
1.04 


3 50 + °-^6 

O.OU _ 0.98 


xVndf 


233/219 


232/219 



AG < 




Prog. Ref. [27] 


Prog. Ref. [28] 




0.320 ± 0.009 


0.328 ± 0.009 




1 00+ 15 

1.381 oil 


1 rio + 13 


/5s 


A nn+ 29 

4.081 0:27 


, ^_ + 25 

4.05 1 ^oil 


7s 






Vg 


-0.3091 H44 


-0 192 + o-o^"^ 
u.±yz _ 0.109 




0.391 ^8 


n 9q +0.063 
- 0.47 


/3g 


13.91 li 


13.81^1 


"3 


-0.212 ±0.027 


-0.209 ±0.027 


/93 


9 44+ 0.11 
0.10 


2.40 




43 + 

'J-^'J - 0.16 


0.383 1 




f.A+ 0.55 
"J-^^- 0.54 


q qo + 0.33 
o.oy _ 0.39 


xVndf 


247/219 


247/219 







as 




7E 


VG 




/3g 


0:3 


P3 


as 








0.581 


0.143 


-0.432 


-0.548 


0.549 




-0.075 


-0.118 


0.030 


-0.008 


«E 


-0.492 




0.648 


0.272 


-0.434 


0.452 




0.053 


0.066 


-0.121 


-0.047 


/3e 


-0.388 


0.877 




0.304 


-0.011 


0.022 




-0.010 


-0.037 


-0.420 


-0.499 


7s 










0.272 


-0.248 




0.088 


0.142 


-0.361 


-0.025 




0.277 


-0.221 


-0.130 






-0.978 




0.082 


0.066 


0.071 


0.067 




0.162 


-0.052 


0.012 




0.835 






-0.087 


-0.070 


-0.069 


-0.063 


/3g 


0.148 


-0.039 


0.025 




0.814 


0.935 












03 


-0.012 


0.008 


-0.032 




0.078 


0.006 


0.053 




0.788 


-0.023 


-0.020 




-0.104 


0.067 


0.037 




0.060 


0.003 


0.023 


0.793 




-0.017 


-0.013 


as 


-0.105 


-0.175 


-0.276 




0.171 


0.099 


0.219 


-0.036 


-0.016 




0.832 




-0.137 


0.033 


-0.211 




0.118 


0.063 


0.138 


-0.044 


-0.026 


0.821 





Table 3: Top: Values of the parameters obtained from the QCD analysis at = 3 (GeV/c)^ in fits with 
AG > and AG < with the two programs. The quoted errors correspond to one a and have been 
obtained from the MINOS analysis [34]. The strongly asymmetric errors obtained for some parameters 
are due to the positivity constraints applied in the fits. Bottom: Correlation matrices for the fits by the 
program of Ref. [27]. The triangles above and below the diagonal correspond to the fits with AG > 
and AG < 0, respectively. The "-" symbols correspond to parameters which are fixed in one of the fits. 
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Range in x 


COMPASS data evolved to Q"^ = 3{GeY/cf using 


fits of 


COMPASS fits (prog. [27]) 


BB [16] 


LSS [18] 


AG > 


AG < 


[0.004, 0.7] 


0.0455 


0.0469 


0.0469 


0.0511 


[0.7, 1] 


0.0014 


0.0008 


0.0011 


0.0010 


[0, 0.004] 


-0.0040 


-0.0029 


-0.0014 


0.0004 


[0, 1] 


0.0430 


0.0448 


0.0466 


0.0525 



Table 4: Contributions to Tf (Q^ = 3(GeV/c)2) from different kinematic regions. The values in the 
first line are the COMPASS results evolved according to different fits and integrated over the measured 
X range. The second and third lines show the corresponding high and low x extrapolations. 
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Figure 1 : Fraction of inclusive, semi-inclusive, and 
calorimetric triggers as a function of x. Events are 
counted with the weight they carry in the asymme- 
try calculation. 



0.001 




Figure 2: Difference between asymmetries for in- 
clusive and hadronic Monte Carlo events in the 
kinematic range covered by the purely calorimetric 
trigger. 
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Figure 3: The asymmetry Af{x) as measured in COMPASS and previous results from SMC [5], HER- 
MES [4], SLAC E143 [2] and E155 [3] at > l(GeV/c)2. The SLAC values of gi/Fi have been 
converted to Ai and the E155 data corresponding to the same x have been averaged over Q^. Only statis- 
tical errors are shown with the data points. The shaded areas show the size of the COMPASS systematic 
errors. 




Figure 4: Values of xgf{x). The COMPASS points are given at the (Q^) where they were mea- 
sured. The SMC points have been moved to the of the corresponding COMPASS points. 
Only statistical errors are shown with the data points. The shaded band at the bottom shows 
the COMPASS systematic error. The curves show the results of QCD fits with AG > and 
AG< 0. 
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X 

Figure 5: The COMPASS values of evolved to Q"^ = 3(GeV/c)^. The open triangles at low x corre- 
spond to Q"^ > 0.7(GeV/c)2, the other symbols to > l(GeV/c)2. Results of QCD fits are shown by 
curves. In addition to our fits (AG > and AG < 0) the curve obtained with three published polarised 
PDF parameterizations (Bliimlein and Bottcher, GRSV and LSS05) [15] is shown. These parameteriza- 
tions lead almost to the same values of (x, = 3(GeV/c)^) and have been averaged. For clarity 
the data points evolved with different fits are shifted in x with respect to each other Only statistical errors 
are shown. 
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Figure 6: Values of Af as a function of in intervals of x. The solid lines show the results of fits to a 
constant. 



14 




Figure 7: Gluon distribution xAG(x) corresponding to the fits with AG > (left) and AG < (right) 
obtained with the program of Ref. [27]. The dashed, solid and dotted lines correspond to = 1.5, 
3 and 10 (GeV/c)^, respectively. The unpolarised distributions ±xG(x) which were used in the fit as 
constrains for the polarised ones are shown for = 1.5 and 3 (GeV/c)^. 




Figure 8: Strange quark distribution xAs(x) corresponding to the fits with AG > (left) and AG < 
(right) obtained with the program of Ref. [27]. The dashed, solid and dotted lines correspond to = 1.5, 
3 and 10 (GeV/c)^, respectively. The unpolarised distributions ±xs(x) are shown for = 1.5 and 
3 (GeV/c)2. 
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10'^ lO"' 



Figure 9: Distribution of the gluon polarisation AG{x)/G{x) at = 3(GeV/c)^ for the fits with 
AG > and AG < obtained with the program of Ref. [27]. The error bands correspond to the 
statistical error on AG{x) at a given x. The unpolarised gluon distribution is taken from the MRST 
parametrisation [26]. The three data points show the measured values from SMC [31], HERMES [32] 
and COMPASS [33]. Two error bars are associated to each data point, one corresponding to the statistical 
precision and the other one to the statistical and systematic errors added in quadrature. The horizontal 
bar on each point shows the x-range of measurement. 
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